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FOREWORD 


This  Interim  Technical  Report  No.  2  was  prepared  for  presentation 
at  the  Conference  on  the  Prospects  of  Advanced  Fracture  Mechanics  which 
is  to  be  held  in  Delft,  Netherland,  June  24  -  28,  1974,  under  the  joint 
sponsorship  of  the  Metal  Research  Institute  T.  N.  0.  and  Delft  Univer¬ 
sity  of  Technology  in  close  cooperation  with  Lehigh  University. 

The  figures  in  this  report  were  reduced  in  compliance  with  the 
format  specified  by  the  Seminar  organisers. 

Full-size,  8  1/2"  x  11",  figures  can  be  obtained  for  numerical 
computation  by  writing  to  the  undersigned: 


Professor  Albert  S.  Kobayashi 
Department  of  Mechanical  Engineering 
University  cf  Washington 
Seattle,  Washington  98195 


THE  FINDINGS  OF  THIS  REPORT  ARE  NOT  TO  BE 
CONSTRUED  AS  AN  OFFICIAL  DEPARTMENT  OF  THE 
ARM*  POSITION,  UNLESS  SO  DESIGNATED  BY  OTHER 
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STRESS  INTENSITY  FACTORS  FOR  ELLIPTICAL  CRACKS 

A.S.  Kobayashi,  A.N.  Enetanya,  and  R.C.  Shah 

University  of  Washington,  Seattle,  Washington  98195  and 
The  Boeing  Aerospace  Company,  Seattle,  Washington  98124 


ABSTRACT 

The  combined  use  of  a  solution  for  an  elliptical  crack  sub¬ 
jected  to  internal  pressure,  which  is  represented  by  a  polynomial 
expression,  together  with  a  free  surface  solution  for  the  alter¬ 
nating  method  is  reviewed-andmathods  for  improving- the- efficiency 
of  this  numerical  technique  are  proposed.  Numerical  differentia¬ 
tion  for  avoiding  the  derivation  and  computer  programming  of  com¬ 
plex  mathematical  expressions  is  discussed.  A  criterion  for  max¬ 
imum  grid  size  for  the  free  surface  solution  ie  proposed.  Also 
optimum  fittihgs  of  the  firr.t»  second  and  third  order  polynomial 
distributions  of  internal  pressure  are  considered. 


INTRODUCTION 

Post  mortem  failure  analysis  of  fractured  structural  compo¬ 
nents  requires  modeling  of  an  embedded  or  surface  flaw  from  which 
fracture  initiated.  More  often,  these  flaws  can  be  approximated 
by  an  ellipse  and  are  located  in  regions  of  stress  concentrations 
or  thin  cross-section,  such  as  those  shown  in  Fig.  1,  where  the 
stress  gradient  as  well  as  the  complex  geometry  of  the  structure 
cannot  be  Ignored.  The  elliptical  crack  and  its  interaction  with 
its  adjacent  or  intersecting  boundaries  have  thus  been  studied  by 
various  investigators  during  the  past  ten  years. 

The  earliest  paper  on^stress  intensity  factor  of  an  ellipti¬ 
cal  crack  is  that  by  Irwin  who  then  derived  the  approximate  ex¬ 
pression  for  the  stress  intensity  factor  at  the  point  of  deepest 


DEF eCT  IN  A  FLANGE  -  WEB  AT  A  HOLE 

CONNECTION 


(c)  DEEP  PART- ELLIPTICAL 
SURFACE  FLAW 


id)  EMBEDDED  ELLIPTICAL  FLAW 
IN  A  THIN-WALLED  PRESSURE 
VESSEL 


FIGURE  I  TYPICAL  ELLIPTICAL  FLAWS  IN  STRUCTURAL  COMPONENTS 


penetration  of  a  shallow*  semi-elliptical  crack  in  a  tension  plate. 
Among  the  several  solutions  on  the  "surface  flaw  problem"  as  well 
as  embedded  flaws  which  have  been  published  to  date,  the  mathemat¬ 


ically  more  rigorous  solutions  are  all  based  g11  the  alternating 
technique  suggested  by  Kantorovich  and  Krylov  .  This  procedure 
-was- first  applied  to_ crack. .problems  by  SmithJ  who  determined  the 
stress  intensity  factors  of  a  semi-circular  crack  in  a  thick  beam 


in  tension,  bending,  and  subjected  to  rapid  heating  or  cooling. 


This  method  of  approach  was  then  used  by  Smith  and  his  coworkers  to 
analyze  a  family-of  problems  involving  circular  flaws  and  part- 
circular  t'laws^’  *  .  Recently  Hartranft  and  Sin  reexamined  one 
of  Smith's  original  solutions with  an  added  singularity  func¬ 


tion  which  increased  the  accuracy  of  their  solution. 


The  disadvantage  of  using  the  circular  crack  solution  is  its 
geometry  which  cannot  be  used  to  model  oblong  elliptical  or  semi¬ 
cracks  such  ns  those  shown  in  Figs.  In  and  ld£*.  For  better  mod¬ 
eling  of  elliptical^laws,  Shah  and  Kobayashi  used  the  potential 
function  of  Segedinx  to  derive  a  polynomial  representation  of  the 
internal  pressure  on  an  elliptical  crack.  This  solution  was  then 


used  in  the  alternating  method  to  determine  the  stress  intensij^ 
factors  of  embedded  elliptical  flaws  under  nonuniform  loadings 


Shallow  crack  means  that  the  crack  depth  is  less  than  half 
of  the  plate  thickness. 

Part  elliptical  cracks,  such  as  those  shown  in  Figs,  lb  and 
lc,  can  be  approximated  by  part  circular  crack  as  shown  in 
Ref.  8.  Part  elliptical  crack  with  large  aspect  ratio,  i.e., 
b/a  <_  0.2,  cannot  be  conveniently  approximated  by  a  part 
circular  crack. 
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as  well'  as  to  estimate  the  stress  intensity  factors  fcr  semi-  _ . 

elliptij^l  surface  flaws1  .  Prior  to  the  above  development,  Kassir 
and  Sili  determined  the  9tress  intensity  factor  of  an  elliptical 
crack  under  linearly  varying  internal  pressure  and  uniform  shear 
tractions.  Smith3-  extended  Kassir  and  Sih  solution-1  following 
the  procedure  described  in  Ref.  9  to  derive  a  polynomial  distribu¬ 
tion  of  shear  factions  prescribed  on  an  elliptical  crack.  In 
addition,  Smith3-  has  recently  used  this  procedure  to  determine  the 
stress  intensity  factor  of  a  surface  flaw  in  a  thin  tension  plate. 

The  solution  procedure  involving  an  elliptical  crack  is  con¬ 
venient  for  modeling  actual  flaw  geometries.  The  small  number  of 
derivable  terms  in  the  polynomial  representation  of  the  internal 
pressure,  however,  limits  the  accuracy  in  which  a  complex  residual 
pressure  distribution  on  the  crack  surface  can  be  fitted.  Unfor¬ 
tunately  numerical  convergence  of  the  alternating  method  is  sensi¬ 
tive  to  this  fitting  and  thus  the  accuracy  of  the  entire  solution 
procedures  hinges  on  the  accuracy  of  fitting.  In  particularly,  nu¬ 
merical  difficulties  in  fitting  residual  stress  are  encountered  in 
surface  flaw  problems  whera  a  portion  of  the  ellipse  protrudes 
through  the  free  surface. 

The  purpose  of  this  paper  is  to  discuss  numerical  problems 
associated  with  the  fitting  process  and  to  suggest  procedures  which 
reduce  some  of  the  immense  amount  of  theoretical  derivations  and 
computer  programming  associated  with  this  numerical  analysis. — — 


THEORETICAL  BACKGROUND 

Two  solutions  are  needed  in  applying  the  alternating  method  to 
fracture  mechanics  problems  Involving  elliptical  cracks.  One  solu¬ 
tion  involves  a  semi-infinite  solid  and  the  other  solution  involves 
an  elliptical  crack  both  subjected  to  variable  surface  tractions. 

For  eJ.Jjijgf|ical  crack  problems  involving  flat  surfaces,  Love's  solu¬ 
tion  ’  is  used  for  the  free-surface  solution  of  the  former. 
Nu^ejjigaJ^pjgcedures  related  to  Love's  solution  are  well  document¬ 
ed-3’  *  ’  ’  and  will  not  be  repeated  here.  A  discussion,  however, 

on  the  tolerable  maximum  size  of  a  rectangle  for  erasing  the  resid¬ 
ual  surface  traction  on  the  flat  surface  will  be  presented  later. 

In  the  following,  a  brief  review  of  the  basic-equations  for  an 
elliptical  crack  with  prescribed  internal  pressure  ,  which  can  be 
represented  by  a  polynomial  in  terms  of  Segedin's  potential  func¬ 
tion3-  ,  is  presented. 

Consider  an  elastic  solid  containing  an  elliptical  crack,  as 
shown  in  Fig.  2,  which  is  located  in  the  plane  2»0,  and  is  opened 
by  applying  an  internal  pressure,  p(x,y),  symmetrically  to  both  sur¬ 
faces  of  the  crack.  The  boundary  conditions  for  this  problem  are: 
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FIGURE  2.  ELLIPTICAL  CRACK  IN  A  FINITE  SOLID 
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For  the  vanishing  shear  stresses  on  z=0  plane  and  in  the  ab¬ 
sence  of  body  forces,  Navier's  equations  of  equilibrium  are  satis¬ 
fied  by  the  harmonic  function  of  $.  The  displacement  perpendicular 
to  the  crack  plane  in  terms  of  <f>  is, 


3z2 


2  (1  -  n) 


it 

3z 


(2) 


where  n  is  the  Poisson's  ratio.  The  stress  components  necessary 
for  solving  the  surface  flaw  problem  are  given  by, 
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where  G  is  the  shear  modulus. 


(3b) 

(3c) 

(3d) 

(3e) 


Due  to  practical  limitations  in  deriving  additional  terms  in 
the  polynomial  representation  of  the  pressure  distribution,  p(x,y), 
only  terms  up  through  the  third  powers  in  x  and  y  were  derived  as, 
3  i  i 

p(x,y)  -  l  A.,  x  yJ  .  (4) 

i, j“Q  J 

The  harmonic  function  which  satisfies  the  necessary  boundary  condi¬ 
tions  was  derived  from  Segedin's  potential  function.  Details  of 
the  harmonic  function  and  its  derivatives  are  described  in  Refs.  9 
and  14. 


Finally,  the  stress  intensity  factor  associated  with  4  can  be 
obtained  through  a  procedure  described  by  Irwin2®  or  Kassir  and 
Sih  .  Again  details  of  this  procedure  are  described  in  Refs.  9 
and  14 . 


Finite  Difference  Calculation  of  Stresses 

As  mentioned  earlier,  the  derivable  terms  in  the  Eq.  4  were 
limited  to  a  third  order  polynomial  in  x  and  y  due  to  practical 
considerations.  Some  of  this  difficulty  can  be  removed  if  the  par¬ 
tial  derivatives  in  Eq.  3  can  be  replaced  by  numerical  differenti¬ 
ations.  Analytical  expressions  of  32<j>/3x2  etc,,  which  are  given  in 
the  Appendix  of  Ref.  9, were  used  to  obtain  the  third  partial 
derivatives  of  $  with  respect  to  x,  y  and  z.  Analytical  expres¬ 
sions  of  32$/3x2  etc.  were  evaluated  numerically  at  two  locations 
of  x  and  x  +  0.001a  etc.  to  compute  33$/3x3  etc.  These  numerical 
derivatives  were  then  combined  as  per  Eq.  3  to  determine  the 
stresses  of  0^ ,  t  and  x  at  various  values  of  x  and  z  on  vari- 
ious  y  =  constant  planes.  These  numerical  results  were  found  to 
agree  within  three  significant  figures  with  corresponding  stresses 
obtained  previously3’". 

The  above  numerical  differentiation  procedure  has  been  extend¬ 
ed  to  compute  the  same  third  partial  derivatives  from  analytical 
expressions  of  3<p/ etc.  Preliminary  results  show  that  the  same 
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excellent  agreement  between  stresses,  which  were  computed  numeri¬ 
cally  atul  -by  -  the  closed  form  procedure,  exists  except  for  regions 
of  high  stress  gradients  where  the  elliptical  crack  front  inter¬ 
sects  with  the  y  «  const,  plane,  e.g.  y**G  plane,  in  Fig.  2.  Effort 
is  currently  underway  to  determine  whether  the  average  surface 
traction  acting  on  the  small  region  surrounding  such  stress  singu¬ 
larity  can  exhibit  numerically  the  same  effect  as  the  residua] 
singular  surface  traction  as  far  as  the  free  surface  solution  is 
concerned.  If  such  replacement  is  possible,  then  the  average 
stresses  determined  by  the  second  numerical  differentiations  of 
3<)73x,  etc.  can  also  be  used  to  compute  the  residual  surface  trac¬ 
tions  on  the  bounding  free  surfaces. 


Maximum  Permissible  Rectangle  Size 


the  second  step  of  the  alternating  technique  is  to  eliminate 
the  residual  surface  traction  on  the  bounding  free  surfaces.  A 
common  procedure  is  to  use  Love's  solution^8*19  for  a  half  space 
with  prescribed  uniform  surface  tractions  on  a  rectangle  in  the 
bounding  plane.  When  the  alternating  technique  was  first  used  to 


ago,  the  total  number 
Through  numerical 


analyze  surface  flaw  problems  several  year§ 
of  such  rectangles  was  of  the  order  of  540- 

experimentation,  Smith  has  gradually  reduced  the  numbers  of  these 
rectangles  to  180°  and  then  to  711  .  Hartranft  and  Sih,  on  the 
other  hand,  increased  the  number  of  rectangles  to  1744  which  to -r 
gether  with  the  increased  number  of  Fourier  series  used,  resulted 
in  computing  time  of  8  hours  to  solve  one  problem.  The  purpose  of 
this  section  is  to  determine  the  maximum  permissible  rectangle  size 
which  will  yield  reasonably  accurate  results  of  this  second  steo  of 
of  the  alternating  procedure. 


Since  the  continuously  varying  residual  surface  tractions  on 
the  bounding  plane  are  to  be  replaced  by  uniform  surface  tractions 
over  each  of  the  rectangles  forming  the  bounding  plane,  the  prob¬ 
lem  reduces  to  the  difference  in  the  resultant  stresses  due  to 
prescribed  uniform  or  variable  surface  tractions  on  the  rectangle. 
This  difference  in  resultant  stresses  can  be  estimated  by  examining 
the  stresses  in  a  half  space  due  to  (1)  a  linearly  varying  stress 
distribution  over  the  rectangle  on  the  bounding  plane,  and  (2)  a 
uniform  stress  distribution  over  the  rectangle.  The  uniform  stress 
distribution  is  in  equilibrium  with  the  linearly  varying  stress 
distribution  as  shown  in  Fig.  3.  The  stress  distributions  at  any 
point  in  a  semi-infinite  plate  due  to  a  linearly  varying  or  a  uni¬ 
form  surface  traction,  Oyy  ,  over  a  segment  2d  as  shoyji  in  the 
legend  of  Fig.  3  are  given  by  the  following  equations 


a  ** 
zz 


1°  (£  ltl  xltlgri),.2  +  .Jr.teid).,  z+d  -1  ,._2dy - } 

7T  ld  y2+(z+d)2  T  y2+(Z-d)2  d  y2+22-d2; 

(5a) 
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FIGURES.  STRESS  DISTRIBUTION  AT  Z  *  d  IN  A  HALF  PLANE  DUE 
TO  UNIFORM  AND  LINEARLY  VARYING  PRESSURES  ACTiNG 
ON  THE  BOUNDING  LINE. 


T 

zy 


_  - z  tan-l  __2dy - 

71  y2+(z-d)2  ^  y2+z2-d2 


(5b) 


For  a  uniform  surface  traction  prescribed  over  a  segment  of 
these  stresses  become^1- 

0  -  -  ^  (tan"1  — - } 

zZ  1T  y2+z2-d2  iy2+(z+d)2]  [y2+(z-d) 2] 

£o  _ _ 4dyz _ 

Tzy  ’  "  *  [y2+(z-i-d)2][y2+(z-d)2]  ' 


2d, 

(6a) 

(6b) 


Figure  3  shows  the  distributions  of  o2Z  and  Tzy  at  distances 
of  z  «  d  from  the  center  of  the  load  segment  for  linearly  varying 
and  uniform  normal  tractions  acting  on  the  segment.  Although  the 
stresses  due  to  the  two  loads  differ  by  a  factor  of  2  at  y  =  0  in 
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Pig.  3,  this  difference  when' superimposed  on  to  ocher  stresses  due  co 
the  common  uniform  load  acting  on  the  rectangle,  reduces  the  over¬ 
all  differences  to  less  than  5%  in  actual  computations.  Further 
evaluations  of  Eqs.  5  and  6,  for  most  parts  of  the  elliptical 
cr.  k,  show  that  the  differences  in  both  normal  and  shear  stresses 
on  the  crack  plane  due  to  the  above  two  surface  tractions  rapidly 
diminish  at  a  distance  which  is  half-wid^h  away  from  the  load 
segment. 

The  above  results  indicate  thar.  the  size  of  the  rectangle  can 
be  equal  to  its  closest  distance  to  the  crack  plane  or  to  the  other 
bounding  surfaces  for  the  case  of  a  finite  thickness  solid,  which¬ 
ever  is  smaller.  Similar  criterion  was  found  to  hold  for  ♦•angen- 
tial  surface  traction  acting  on  the  bounding  free  surface.  The 
number  of  necessary  rectangles  can  be  reduced  systematically  fol¬ 
lowing  this  criterion  thus  resulting  in  substantial  reduction  in 
computing  time. 


Residual  Pressure  Distribution  on  Semi-Elliptical  Crack 

Another  recognized  drawback  of  the  alternating  technique  based 
on  the  elliptical  crack  potential  function  is  in  its  inability  to 
approximate  the  residual  pressure  distribution  on  the  crack  surface 
with  the  limited  third  order  polynomial.  This  drawback  is  for  all 
practical  purposes  nonexistent  for  totally  embedded  elliptical 
flaw,  such  as  those  shown  in  Figs,  la  and  Id,  since  the  residual 
pressure  distribution  varies  smoothly,  A  third  degree  polynomial 
can  be  fitted  within  5  percent  of  such  residual  pressure  distribu¬ 
tion  and  thus  the  numerical  procedure  converges  rapidly  with  typi¬ 
cal  incremental  change  in  the  stress  intensity  factor  to  be  less 
than  0.2  percent  of  the  original  stress  intensity  factor  after  2  to 
5  iterations. 

For  a  surface  flaw  where  the  elliptical  circle  penetrates  the 
free  surface,  this  fitting  problem  becomes  acute  since  the  residu¬ 
al  pressure  on  the  crack  surface  normally  reaches  a  maximum  value 
in  the  vicinity  of  the  free  bounding  surface.  The  physical  crack 
obviously  does  not  continue  beyond  the  bounding  surface  while  the 
elliptical  crack,  used  in  the  alternating  method,  does  protrude  in¬ 
to  this  free  space.  In  order  to  represent  the  free  bounding  sur¬ 
face,  Hartranft  and  oltr  have  suggested  the  use  of  only  even  terms 
of  polynomial  which  reduces  Eq.  4  to  azz  *  £qq  +  A2q  x2  +  Aq2  y2  to 
which  the  residual  pressure  distribution  must  be  fitted.  The  re¬ 
sult  is  an  incomplete  erasure  of  the  residual  pressure  distribu¬ 
tion  resulting  in  slow  convergence  or  possible  divergence  of  the 
iteration  process  in  the  alternating  procedure. 

2 

In  a  recent  paper,  one  of  the  authors  showed  that  the  stress 
intensity  factor  for  a  circular  crack  of  radius  b  embedded  in  an 
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infinite  solid  and  subjected  to  linearly  varying  pressure  distribu¬ 
tion  of  o.£Z  »  p0(l  -  ^-)  was  approximately  5  percent  higher  than 
that  of  a  semi-circular  surface  crack  subjected  to  the  same  pres¬ 
sure  distribution  between  0  '  y  <  b.  This  empirical  finding  indi¬ 
cated  that  the  stress  intensity  magnification  factor  due  to  the 
free  bounding  surface  had  an  effect  similar  to  continuing  the  lin¬ 
early  varying  pressure  distribution  in  the-b  £  y  0  region  of  the 
fictitious  portion  of  the  semi-circular  crack.  This  approximation 
would  obviously  result  in  a  significant  error  in  estimation  of  the 
stress  Intensity  magnification  factor  If  a  similar  procedure  were 
used  to  continue  into  the  region  of  -b  <,  y  <,  0,  the  nonlinear  pres¬ 
sure  distribution  acting  on  a  semi-circular  crack. 

The  above  empirical  finding  led  to  further  numerical  experi¬ 
mentation  for  arriving  at  an  optimum  pressure  distribution  on  the 
crack  surface  extended  into  the  region  of  -b  <_  y  0  of  the  ellip¬ 
tical  crack.  Three  pressure  distributions  of  linear,  quadratic, 
and  cubic  with  the  maximum  pressure  at  the  free  surface,  as  shown 
in  Fig.  4  were  considered.  For  twc -dimensional  problems  the  stress 
intensity  factors  associated  with  an  edge  crack  with  linear,  quad¬ 
ratic,  and  cubic  pressure  distributions,  shown  in  Figs.  4a,  4c,  and 
4e,  can  be  computed  by  using  the  results  of  Stallybrass2^,  In  order 
to  obtain  the  same  stress  intensity  factor  at  the  crack  tip  of  a 
totally  embedded  through  crack,  various  linear  distributions  of 
pressures  were  prescribed  on  the  other  half  of  the  embedded  crack, 
i.e.,  -b  <,  y  <,  0.  It  was  found  that  preoaure  distributions  shown 
In  Figs.  4b,  4d,  and  4f  yielded  stress  intensity  factors  which  were 
reasonably  close  to  those  for  the  edge  crack*.  These  pressure  dis¬ 
tributions  were  then  prescribed  onto  the  other  half  of  the  ell  ',ptl- 
cal  crack,  i.e.,  -b  <  y  <  0,  where  half  of  the  crack  was  loaded  by 
linear,  quadratic,  or  cubic  pressure  distributions.  Six  terms  in 
Eq.  4  or  ozg  “  Aqq  +  Ag^y  +  A2QX2  +  Agjy2  +  A2ix2y  +  Ag3y3  were 
then  least  square  fitted  to  this  prescribed  pressure  distribution 
on  the  crack  surface. 

The  above  procedure  of  least  square  fitting  was  also  applied 
to  the  two-dimensional  problem  illustrated  in  Fig.  4.  For  this 
purpose,  %  pressure  of  azz  •  Aq  +  A^y  +  A£y2  +  A3y3  was  fitted  to 
the  3  pressure  distributions  illustrated  in  Figs.  4b,  4d,  and  4f. 
Table  I  shows  the  fitted  stress  distributions,  stress  Intensity 
factors  at  the  right  crack  tip,  and  the  corresponding  stress  in¬ 
tensity  factors  for  edge  cracks. 


* 


Stress  intensity  factors  at  the  right  crack  tip  were  computed 
by  the  procedure  described  in  Ref.  24. 
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K  *  0,4385 Po'/ifb  K  ■*  0.4588  p  </Vb 


(a)  (b) 


K  «0.3»59|V^£>  K  *0.3446PoV^l 

(c)  (d) 


K. 0.6071  Po^rB  K  *  0.6808  p0/W 

(«)  (f) 


PIGURE  4  STRESS  INTENSITY  FACTOR  EOS  AN  EDGE 
CRACK  SUBJECTED  TO  VARIABLE  PRESSURE 
LOADING 


Differences  between  the 
respective  stress  in¬ 
tensity  factors  can 
obviously  be  reduced  by 
further  numerical  exper¬ 
imentation  with  various 
pressure  distributions 
prescribed  in  the  region 
of  -b  <  y  <_  0.  The  pur¬ 
pose  of  thTs  study,  how¬ 
ever,  was  to  demonstrate 
that  reasonable  results 
can  be  obtained  for  the 
three  cases  of  prescribed 
linear,  quadratic,  and 
cubic  pressure  distribu¬ 
tions  with  the  same 
additional  pressure  in 
the  region  of  -b  _<  y  <,  0 
and  thus  no  further 
attempt  was  made  to 
improve  the  accuracy  of 
the  resultant  stress 
intensity  factors. 


TABLE  I 


Prescribed 

Pressure 

Fitted  Pressure 

Kj (Right 
crack  tip) 

Kj  (Edge  Crack) 

Linearly 

Varying 

Fig.  4b 

o  -0. 9960-0. 8714y 

28  -0 . 4480y ^+0 . 3350y 3 

0.4636 

P  M 

*Q 

0.4385  pQYnb 

Quadratically 

Varying 

Fig.  4d 

o  -1. 0042-0. 1843y 

2  -0.4937y2-0,4094y:' 

0.5169 

p  /rfb 

*0 

0.5950  pQv^b 

Cubicully 

Varying 

Fig.  4f 

o  “0.9965+0. 1508y 

22  -0.4480y2-0.7276y' 

0.5808 

P  Ab 
*  0 

0.6071  pQnib 

Figure  5  shows  the  normalised  stress  intensity  factor  of  a 
circular  crack  for  the  pressure  distribution  shown  in  Fig.  4b. 
Poisson's  ratio  for  this  analysis  was  set  to  n  *  0.3.  The  poly¬ 
nomial  pressure  distribution  fitted  within  2%  of  the  prescribed 
bilinear  pressure  distribution.  As  shown  in  Fig.  5*  the  stress 
intensity  factor  thus  obtained  differed  as  much  as  9%  from  the 
corresponding  stress  intensity  factor  for  a  semi-circular  surface 
crack  with  a  prescribed  linearly  varying  load.  Part  of  this  dif¬ 
ference,  i.e.,  approximately  2%,  can  be  attributed  to  the  differ¬ 
ences  in  the  Poisson's  ratios  used  in  the  two  analyses®. 

It  is  interesting  to  note  that  the  residual  normal  surface 
tractions,  Oyy,  on  the  y  *  0  plane  for  the  problem  illustrated  in 
the  legend  of  Fig.  4  are  nearly  equal  to  the  corresponding  resid¬ 
ual  normal  tractions  for  a  circular  crack  subjected  to  uniform 
pressure  of  p0.  It  appears  then  that  the  additional  residual  tan¬ 
gential  surface  tractions,  Ty23y»o»  Problem  act  to  counter¬ 

act  the  crack  opening  moments  of  the  residual  normal  surface  trac¬ 
tions  thus  causing  the  resultant  stress  intensity  factor  to 
approach  that  of  a  semi-circular  crack  in  a  semi-infinite  solid  and 
subjected  to  linearly  varying  pressure. 

The  close  agreement  between  the  stress  intensity  factors  for 
the  two  extremes  of  a  semi-circular  crack  and  a  two-dimensional  edge 
crack  suggested  that  this  procedure  of  prescribing  appropriate 


riGUSE  S.  STRESS  INTENSITY  FACTORS  FOR  A  CIRCULAR  CRACK  SUBJECTED  TO  A  COMBINED  TWO 
LINEAR  VARYING  PRESSURE  DISTRIBUTIONS  ANOA  semi. circular  SURFACE  CRACK  WITH 
linearly  VARING  load 
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pressure  distribution  in  the  region  of  -b  <_  y  £  0  of  an  elliptical 
crack  can  be  used  to  estimate  the  stress  intensity  factors  for 
semi-elliptical  surface  flaws  subjected  to  linearly  varying  pres¬ 
sure  distribution.  The  estimated  stress  intensity  magnification 
factors  thus  obtained  along  the  periphery  of  a  semi-elliptical 
crack  subjected  to  linearly  varying  pressure  distribution  are 
shown  in  Fig.  6.  The  stress  intensity  factor  in  the  vicinity  of 
6*0,  i.e.,  where  the  semi-elliptical  crack  meets  the  free  bound¬ 
ing  surface,  is  not  shown  since  the  stress  singularity  at  this 
point  is  believed  to  vanish  following  the  analysis  by  Hartranft 
and  Sih7.  The  error  bounds  of  this  estimation  are  indicated  in 
Fig.  6  by  the  two  known  stress  intensity  factors  for  b/a  *  0  and 


1.0.  The  overestimation  in  the  stress  intensity  magnification  fac¬ 
tors  at  the  deepest  penetration  of  the  ellipse  or  8  »  90°  is  con¬ 
sistent  with  the  corresponding  two-dimensional  problem  shown  in 
Table  I.  For  a  better  estimate  of  the  stress  intensity  magnifica¬ 
tion  factor  at  this  position,  all  values  could  be  lowered  7%  fol¬ 
lowing  the  results  In  Table  1. 

The  same  procedure  was  used  to  estimate  the  stress  intensity 
magnification  factors  for  semi-elliptical  cracks  subjected  to 
quadratic  and  cubic  distributions  of  pressure.  Again,  the  fitted 
polynomial  distribution  of  pressure  fitted  within  2  percent  of  the 
prescribed  pressure.  Unlike  the  previous  case  of  linearly  varying 
pressure  distribution,  the  stress  intensity  factors  for  quadratic 
and  cubic  pressure  distribution  on  a  semi-circular  surface  flaw 
are  not  known.  The  obtained  stress  intensity  factors,  however, 
were  consistent  with  known  theoretical  solutions  at  one  extreme  of 
b/a  “  0,  i.e. ,  for  an  edge  crack. 


DISCUSSION 

The  stress  intensity  factor  of  a  semi-elliptical  crack  in  a 
semi-infinite  solid  subjected  to  uniaxial  tension  can  be  estimated 
from  the  results  of  Fig.  5.  Figure  7  shows  the  pressure  distribu¬ 
tion  for  erasing-the  residual,  normal  stresses  on_the  free  surface 
at  y«0  of  a  semi-infinite  solid  after  the  first  iteration  in  the 
alternating  method.  An  equivalent  linearly  varying  pressure  dis¬ 
tribution  as  shown  in  Fig,  7,  which  is  nearly  in  equilibrium  with 
the  rapidly  varying  pressure  distribution  is  now  considered.  The 
stress  intensity  factor  for  a  semi-elliptical  crack  subjected  to 
the  equivalent  linearly  varying  pressure  distribution  is  then  esti¬ 
mated  from  Fig.  6.  The  superposition  of  this  stress  intensity 
factor  to  the  original  stress  intensity  factor  of  an  elliptical 
crack  in  an  infinite  solid  and  subjected  to  uniform  pressure  then 
yields  an  estimate  of  the  front  surface  magnification  factor  for  a 
semi-elliptical  crack  in  a  semi-infinite  solid  as  shown  in  Fig.  8. 
It  should  also  be  noted  that,  as  evident  from  the  equivalent  pres¬ 
sure  distributions  In  Fig.  7,  this  procedure  will  yield  underesti¬ 
mated  stress  intensity  magnification  factors  close  to  the  free  sur¬ 
face.  Nevertheless,  for  a  semi-circular  crack,  the  estimated  val¬ 
ue  of  Mp  ■  1.13  at  an  angular  orientation  of  6  ■  20®  is  suffici¬ 
ently  close  to  the  original  results  of  Mp  «  1.1?.  obtained  by 
Smith  et  al. 

The  estimated  front  surface  stress  intensity  magnification 
factors  together  with  those  by  Smith  and  Alavi^ ,  Kobayashi  and 
Moss^S,  and  Nisitani  and  Murakami*®  are  shown  in  Fig.  9.  The  lat¬ 
ter  results  by  Nisitani  and  Murakami,  which  were  computed  for 
Poisson's  ratio  of  n  »  0,  were  corrected  for  n  ■  0.25,  as  shown  in 
Fig.  9  following  an  investigation  on  the  effects  of  Poisson's 
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DISTANCE  ~  OR  ■£ 

FIGURE  7.  PRESCRIBED  PRESSURE  DISTRIBUTION  ON  A  SEMI- 

ELLIPTICAL  CRACK  SURFACE  AFTER  FIRST  ITERATION 
IN  THE  ALTERNATING  METhOO,  POISSON  RATIO  1J-0.3 


ratio  by  and 

Kobayashi“ ' .  Also 
shown  in  Fig.  9  is 
the  single  result  by 
Hartranft  and  Sih 
for  n  «  0.30? . 
Although  the  new 
estimates  of  the 
front  surface  magni¬ 
fication  factor 
appear  to  be  in 
closer  agreement  with 
the  results  by 
Smith  and  Alavi5, 
whan  an  adjustment 
is  made  to  account 
for  the  possible 
overestimation  dis¬ 
cussed  previously i 
the  adjusted  front 
surface  stress  in¬ 
tensity  magnifica¬ 
tion  factor  will  move 
closer  to  the  analyt¬ 
ically  more  rigorous 
results  by  Nisitani 
and  Murakami. 

The  above  men¬ 
tioned  approximate 
procedure  does  not 
eliminate  the  diffi¬ 
culties  in  fitting  a 


polynomial  pressure 

distribution.  Much  of  these  difficulties,  however,  can  be  reduced 
by  a  shrewd  use  of  the  limited  terms  in  the  polynomial  using  nu¬ 
merical  procedures  akin  to  over-and-under  relaxations. 


The  alternating  procedure  with  optimum  rectangle  spacing  on 
the  free  bounding  surface  and  the  pressure  distribution  in  the 
region  of  -b  <_  y  <_  0  of  the  elliptical  crack  are  being  used  to 
solve  practical  problems  involving  semi-elliptical  surface  flaws. 
Hopefully,  these  results  will  be  available  for  presentation  at  the 
Conference  of  Prospects  of  Fracture  Mechanics. 
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